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ABSTRACT 

We compare the results of IGrimm" et all 112003) and iRanalli et all l|200.^ on the L x - 
SFR relation in normal galaxies. Based on the Lx — stellar mass dependence for 
LMXBs, we show, that low SFR ( < 1 MQ/year) galaxies in the Ranalli et al. sample 
are contaminated by the X-ray emission from low mass X-ray binaries, unrelated to 
the current star formation activity. 

The most important conclusion from our comparison is, however, that after the 
data are corrected for the "LMXB contamination" , the two datasets become consis- 
tent with each other, despite of their different content, variability effects, difference in 
the adopted source distances, X-ray flux and star formation rate determination and 
in the cosmological parameters used in interpreting the HDF-N data. They also agree 
well, both in the low and high SFR regimes, with the predicted Lx _ SFR dependence 
derived from the parameters of the "universal" HMXB luminosity function. This en- 
couraging result emphasizes the potential of the X-ray luminosity as an independent 
star formation rate indicator for normal galaxies. 
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1 INTRODUCTION 



Based on Chandra observations of nearby star form- 
ing galaxies and studies of high mass X-ray bina- 
ries (HMXB) population in the Milky Way and SMC, 
IGrimm. Gilfanov fc Sunvaevl (120031) proposed recently that 
high mass X-ray binaries can be used as a star formation 
rate (SFR) indicator. They found, that in a broad range of 
star formation regimes and rates the X-ray luminosity dis- 
tribution of HMXBs can be approximately described by a 
"universal" luminosity function - a power law with the slope 
of ~ 1.6 and a cut-off at lg(Lx) ~ 40.5, which normaliza- 
tion is proportional to the SFR. As the 2-10 keV luminosity 
Lx of a normal galaxy with sufficiently high SFR/M* ratio 
(M 4 - total stellar mass) is dominated by the emission from 
high mass X-ray binaries, the X-ray luminosity can be used 
as a star formation rate indicator for normal galaxies. 

Although the normalization of the luminosity function 
and the number of sources are proportional to the SFR, the 
Lx-SFR dependence is non-linear in the low SFR regime 
and becomes linear only at sufficiently high values of SFR 
(thick solid line in Fig0. This non-linear behavior at low 
SFR values is not related to intrinsic non-linear SFR de- 
pendent effects in the population of the HMXB sources. 
It is rather caused by the fact that the quantity of in- 
terest is a sum of the luminosities of discrete sources - 
Lx,to\ = ^2 k Lx,k, with L x ,k obeying a power law lumi- 
nosity distribution. The non-linear behavior is caused by 



the properties of p(Lx,tot) probability distribution, namely, 
the difference between its expectation mean (average) and 
its mode (mo s t prob able value). This effe ct was discussed in 
IGrimm et all fcOOSt a nd iGilfanovl (I2003T) and will be give n 
a detailed treatment in lGilfanov. Grimm fc Sunvaevl |,2003). 
The position of the break in the Lx-SFR relation is de- 
fined by the parameters of the luminosity function. For par- 
ticular values of the slope and cut-off luminosity found by 
IGrimm et alJ (|2003f), the boundary between non-linear and 
linear regime lies at SFR~ 4.5 M^/year or, equivalently, 
Lx ~ 3 • 10 40 erg/sec. Chandra and ASCA measurements 
of the total X-ray luminosity of a number of nearby star 
forming galaxies were in a good qualitative and quantita- 
tive agreement with the predicted Lx-SFR relation (Fig0 
thick solid curve and filled circles). Moreover, the distant 
star forming galaxies, observed by Chand ra in the Hub- 
ble Deep Field North ferandt et al.l l200ll) at redshifts of 
2 ~ 0.2 — 1.3, also obey the same relation. In the linear high 
SFR regime it is given by: 



SFR[M & /yr] = 



6.7 • 10 39 erg/s 



(1) 



where SFR is the formati on rate of massive stars, M > 
5 Mq . IGrimm et alJ i2003T) pointed out importance of two 
contaminating factors, unrelated to the current star forma- 
tion activity: (i) emission of the central supermassive black 
hole, which even in the low luminosity AGNs can easily out- 
shine X-ray binaries and (ii) contribution of the low mass 
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Fig ure 1. Ly- SFR rela tion. All points from lRanalli et al.l foOOST) 
and lGrimm et all J2003t) . The galaxies with the expected LMXB 
fraction exceeding 50% are plotted as upper limits. The thick 
solid line shows predicted relation between the most probable 
value of Lx and SFR, the shaded area — it's 67% intrinsic spread. 
The straight dashed line shows the expectation mean for Lx, 
which would be obtained if X-ray luminosities of many galaxies 
with similar SFR were averaged. To demonstrate importance of 
the LMXB contribution at low SFR/M* , both HMXB and total 
luminosities are plotted for the Milky Way. 



X-ray binaries, which might be especially important in the 
low SFR regime. 

iRanalli. Comastri fc Setl ]2003l) independently studied 
X-ray luminosity of normal galaxies using the ASCA and 
BeppoSAX archival data and Chandra observations of the 
HDF-N and found a tight correlation between their X-ray, 
radio (1.4 GHz) and FIR fluxes. They suggested that the 2- 
10 keV luminosity of normal galaxies can be used as a SFR 
indicator and derived the relation: 

SFR[M Q /yr] = kcV (2) 

5 • llr y erg/s 

This formula agrees reasonably well with that obtained by 
iGrimm e t alJ (12 003ft for t he high SFR regime, eq.Q. How- 
ever. iRamdliet al.l ]2003T) noted, that the Z/x-SFR relation 
was linear in the entire range of the star formation rates, 
including the low SF R regime, in apparent contradiction to 
IGrimm et ail (120031) results. 

In this L e tter w e compare the lGrimm et al] l|2003h and 
IRanalli et al] |2003) samples of the galaxies. We demon- 
strate, tliat_tlie_2^ay_emission from the low SFR galaxies 
in the IRanalli et all <l2003l) sample is likely to be "contami- 
nated" by low mass X-ray binaries, which are unrelated to 
current star formation activity. After the "LMXB contami- 
nation" is accounted for, the two datasets agree qualitatively 
and quantitatively and are consistent with the Lx~SFR re- 
lation expected on the bas is of the "universal" HMXB lumi- 
nosity function derived bv lGrimni et al] ]2005t . 
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Figure 2. Comparison of the data for local (circles) and HDF- 
N (triangles) galaxi es present both in IGrimm et al] fc003l) and 
IRanalli et alJ l2003t) samples. For each galaxy, it's positions in 
two samples are connected by a broken line with the arrow di- 
rected from G to R. The first segment of each broken line shows 
the effect of the difference in the source distance or cosmological 
parameters, the second segment shows cumulative effect of other 
factors, such as variability and difference in the SFR values. 



2 THE SAMPLES 

In the following we denote IRanalli et al] ((2003) and 
IGrimm et al] ]2003l) samples as R and G correspondingly. 
The data from both samples are plotted together in FigQ 



2.1 The local galaxies 

The two samples, although differently constructed, overlap 
substantially, with 9 galaxies (out of 23 in each sample), 
present in both. The sample R was derived using more rig- 
orously defined construction algorithm. In almost all cases 
the authors adopted diffe r ent d istances and different val- 
ues of SFR. IGrimm et ail <l2003h derived SFR values aver- 
aging the results of several independent estimators based 
on UV, FIR , Hq a nd radio flux measurements, whereas 
IRanalli et al] ((2003) used radio flux measurements. The X- 
ray fluxes were obtained from different observations, some- 
times by different instruments and are, obviously, affected 
by variability of the X-ray emission from the galaxies. For 
some of the g a laxies the X-ray luminosity was calculated by 
I Grimm et al] |2003) as a direct sum of the luminosities of 
compact sources detected by Chandra. 

The Fig|2] compares positions of the galaxies present 
in the both samples in the Lx~SFR plane. Note, that the 
difference in the adopted distances does not have effect at 
high values of SFR where the Lx~SFR relation is linear, but 
it might destroy the correlation in the non-linear low SFR 
regime. 
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2.2 Hubble Deep Field North 

Both iGrimm et all J2003I) and iRanalli et all UpOS) used 
similar selection criteria. Each sample contains seven 
sources, of which six are present in both samples. 
The sources #185 and #148 (according to Table 2 in 
iBrandt et alJ {2001)) are absent from the R and G corre- 
spondingly. The latter was excluded from sample G because 
no 1.4 GHz flux was d etected, with the upper limit of 23 ^Jy 
jRichards et alll998ft . The main difference lies in computing 
the X-ray fluxes and luminosities. IGrimm et al.l i2003T) used 
2-8 keV fluxes from Chandra catalog and K-corrected them 
to 2- 10 keV rest-frame luminosity using the spe ctral indexes 
from lBrandt et alJ feOOlft . IRanalli et al.1 (12003!) derived the 
X-ray count rates in two redshift-corrected energy bands and 
based their final K-correction on the recomputed spectral 
indexes. The following cosmological parameters were used: 
H = 50 km/s/Mpc, q = 0.1 (sample R) and H = 70 
km/s/Mpc, qo = 0.5, A = (sample G). The positions of 
the data points in the Lx-SFR plane are compared in Fig|2] 
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3 LMXB CONTRIBUTION 

Due to long evolutionary time scale, the population of 
low mass X-ray binaries is unrelated to the current star 
formation activity. It is, rath er, proportiona l to the stel- 
lar mass of the host galaxy jGilfanovlEiolh . Hence, the 
X-ray emission from LMXBs can contaminate the Lx- 
SFR relation, as exemplified by the Milky Way galaxy, 
in which the LMXBs contribution exceeds w 90% (FigQ 
IGrimm. Gilfanov fc Sunvaevl l2002ll . Although LMXB and 
HMXB sources can not be easily separated based on the X- 
ray data, and optical identifications are (potentially) avail- 
able only for the most nearby galaxies, the number and 
combined luminosity of LMXBs can be sufficiently accu- 
rately predicted based on the stellar mass of the host galaxy 
jGilfanovi r2003). Thus, relative contributions of LMXB and 
HMXB sources to the X-ray luminosity of the galaxy are 
defined by its position on the SFR— M* plane (Fig|3J. 

Stellar masses of the galaxies were calculated us- 
ing K-band magnitude s from 2MASS Large Galaxy At- 
las llJarrett et aljl200oT) with the color based c orrection to 
the mass-to- light ratio llBell fc de Jond l200lf) . The mass 
of the Milky Way was calculated u sing it s K-band lumi- 
nosity obtained by Malhotra et al ] \ 1996ft from 3D mod- 
eling of the DIRBE data, and assuming the same mass- 
to-light ratio as in M31. The stellar masses of LMC 
and SMC were est imated from their dynamical masses 
jGrimm et alJ [20031. assuming M^ yn /M t — 5. The dis- 
tances for th e galax ies from the sample G are same as in 
IGrimm et alJ J2003) . We re-examined the distance to low 
SFR galaxies. For N GC55 (1.6 Mpc) and M101 (7.2 Mpc) we 
adop ted values from|Puc he 1 Carignan fc Wainscoati \ 1991ft 
and lJurcevic. Pierce fc^agobvl ll2000ft . The distances to 
NGC2403 (3.7 Mpc), NGC2903 (9.5 Mpc), NGC4449 (3.8 
Mpc) and NGC4654 (17.6 Mpc) were est imated from IR 
Tully-Fisher relation IIAaronson et al.ll 982) usin g data from 
iTormen fc BursteirJ i 1995ft and calibration from lSakai et alJ 
i 2000ft . The d istances to other gal axies from R sample are 
the same as in IRanalli et all <2003ft . 

The galaxies from R and G samples are plotted in the 



Figure 3. Location of galaxies from R (open circles) and G (filled 
circles) samples on the SFR— M* plane. The dashed lines corre- 
spond to constant ratio of the the most probable values of HMXB 
and LMXB luminosities estimated from respective average lumi- 
nosity functions with account for non-linear effects of statistics. 
For the points above the solid line, the LMXB contribution ex- 
ceeds 50%. 



SFR— M* plane in Fig|3 along with the contours of con- 
stant Z/hmxb : Ilmxb ratio. The luminosities of LMXB and 
HMXB sources were estimated from their respectiv e aver - 
age luminosity fun ctions obtained by IGrimm et al.l J20 



and lGilfanovl 112003ft . In estimating the LMXB luminosity we 
used average normalization for late type galaxies. Although 
in the limit of large number of sources linear relations hold, 
Llmxb tx M* and Lhmxb oc SFR, the contours are not 
straight lines at M, < 2 ■ 10 10 M w and SFR < 4 M /yr due 
to effects of statistics l|Gilfanov et aljr2003ft . 

As expected, the "LMXB contamination" plays role 
mostly at low SFR values and becomes unimportant at high 
star formation rates. In all but one galaxies from the sam- 
ple R having SFR< 1 MQ/yr, the expected contribution of 
LMXBs exceeds ~ 50%. These galaxies are shown in FigQ 
as upper limits. For two galaxies (NGC55 and NGC2403) 
mostly deviating from the common trend in Fig0 the ex- 
pected LMXB contribution exceeds ~ 70%. 



4 COMBINED SAMPLE AND PREDICTED 
Lx-SFR RELATION 

The Fig0]shows all data from the samples G and R, exclud- 
ing duplicatio ns. The tightnes s of Lx — relation at large 
M, (Fig.14 in lGilfanovlliooil allows one to approximately 
correct observed luminosities for the LMXB contribution. 
This correction was applied for all galaxies, except those 
with M* < 2 • 10 10 M and Lhmxb : Llmxb < 1. The latter 
(NGC55, NGC2403 and IC342) are not plotted in FigH 

The solid curve in Fig|l] shows the predicted Lx~SFR 
relation, calculated using the parameters of the "universal" 
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Figur e 4. Ly— SFR relation - combined data from lRanalli et alJ 
(2003) and|G rimm et alj|2003h . with duplications excluded. The 
Lx values for nearby galaxies were corrected for the LMXB 
contribution estimated from the stellar mass. Three galaxies 
with more than 50% LMXB contribution and small stellar mass, 
M* < 2 ■ 10 10 Mq, for which large intrinsic dispersion of the 
Lx — M* relation precludes accurate estimate of the LMXB lumi- 
nosity, are not plotted. The luminosities for the HDFN and Lynx 
field galaxies were computed for Ho = 70 km/s/Mpc, f2 m = 0.3, 
A = 0.7. The solid and dashed lines and the shaded area are the 
same as in FigUl 

HMXB luminosity function derived bv lGrimm et alJ <l2003t) 
from analysis of five nearby star forming galaxies with best 
known luminosity functions. It corresponds to the mode of 
the probability distribution - the most likely value of the X- 
ray luminosity of a randomly chosen galaxy. The dashed line, 
on the contrary, shows the expectation mean - the value, 
that would result from averaging of the X-ray luminosities 
of many galaxies having similar values of SFR. Due to the 
properties of the probability distribution of the total lumi- 
nosity of a population of discrete sources, Lx,tat = y\ Lx,k, 
these two quantities are not identical in the low SFR limit, 
when the number of sources is small. 

Due to skew ness of the probability distribution 
p{L x ,tot) (Fig.2 in lGilfanovl l2003). large and asymmetric 
dispersion around the solid curve in Figgis expected in the 
non-linear low SFR regime. The probability to find a galaxy 
below the curve is « 12 - 16% at SFR=0.2 - 1.5 M /yr 
and increases to « 30% at SFR=4 — 5 Mq/jx, near the 
break of the Lx-SFR relation. Of course in the linear regime 
(SFR > 10 Mq/yt) it asymptotically approaches ~ 50%, as 
expected. This asymmetry is already seen from the distri- 
bution of the points in Fig0J- at low SFR values there are 
more points above the solid curve, than below. Moreover, the 
low probability high luminosity tail of the p(Lx.tot) distri- 
bution will lead to appearance of galaxies-outliers with sig- 
nificantly larger than expected value of the total luminosity. 
Such galaxies will inevitably appear as the plot is populated 
with more objects. Non-gaussianity of the p(Lx,tot) distribu- 



tion makes least square and \ 2 fitting techniques inadequate 
for analysis of the Lx-SFR relation in the low SFR regime. 



5 CONCLUSION 

We compared res ults of iGrimm et alJ l|2003f) and 
iRanalli et alJ i200ot) on relation of the X-ray luminos- 
ity and the star formation rate in normal galaxies (Fig0 
and El). 

Addressing the discrepancy in th e low SFR regime, w e 
note that six out of seven galaxies from lRanalli et all {2003), 
having SFR < 1 Mq / yr, are likely to be contaminated by the 
X-ray emission from low mass X-ray binaries, having no re- 
lation to the current star formation activity. Furthermore, at 
M* < 10 10 M e and SFR < 1 M /yr, the expected luminosity 
of X-ray binaries does not exceed < 10 39 erg/s. This is com- 
parable or smaller than that of low luminosity AGNs often 
found by Chandra in otherwise apparently normal galaxies. 
The AGN contribution can not be identified and separated, 
unless high angular resolution imaging data are available. 
Secondly, the probability distribution of the total luminos- 
ity of a population of discrete sources, Lx,tot — /, u Lx.fc, 
is significantly non-Gaussian for low values of Lx.tot- This 
should not be ignored when analyzing and interpreting the 
Lx-SFR relation in the low SFR regime. 

The most important conclusion is, however, that af- 
ter the potentially "LMXB contaminated" galaxies are ex- 
cluded, the two datasets become consistent with each other, 
despite of their different content, variability effects, differ- 
ence in the adopted source distances, X-ray flux and star 
formation rate determination and in the cosmological pa- 
rameters used in interpreting the HDF-N data. The ~ 30% 
difference in the calibration of the Lx-SFR relation is in- 
significant considering the number and amplitude of the un- 
certainties involved. They also agree well, both in the low 
and high SFR regimes, with the predicted Lx-SFR depen- 
dence derived from the parameters of the "universal" HMXB 
luminosity function (Fig|3J. This is an encouraging result 
emphasizing the potential of X-ray luminosity as an inde- 
pendent star formation rate indicator. 
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